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Abstract 

We construct a fake exponential Brownian motion, a continuous martingale different from classical 
exponential Brownian motion but with the same marginal distributions, thus extending results of Albin 
and Oleszkiewicz for fake Brownian motions. The ideas extend to other diffusions. 

1 The problem 

Brownian motion is a martingale with centred Gaussian marginals, and variance t. Hamza and Klebaner [5] 
asked, are there other real-valued martingale processes with the same univariate marginals as Brownian 
motion? Since Dupire [5] has shown that Brownian motion is the unique martingale diffusion with 7V(0,f) 
marginals, finding alternative martingales with the same marginals involves relaxing either the continuity 
assumption, or the Markov property. 

Hamza and Klebaner [9] answered their question in the positive, but it turns out that Madan and Yor [17] 
had already provided a construction based on the Azema-Yor solution of the Skorokhod embedding problem. 
Both the Madan and Yor [T7] and Hamza and Klebaner [5] constructions are discontinuous processes. How- 
ever, Albin [I] produced an ingenious solution to the problem with continuous paths, based on products of 
Bessel process. Subsequently, Oleszkiewicz |18j produced a simpler construction which is extremely natural 
and is based on mixing and the Box- Jenkins simulation of a Gaussian variable. Oleszkiewicz also introduced 
the evocative term fake Brownian motion to describe a (continuous) martingale with the same univariate 
marginal distributions as Brownian motion. 

Now that there has been serious study of fake Brownian motion it seems reasonable to ask if there exist 
fake versions of other canonical martingales. In a related context, Hirsch et al [12] give several general 
methods of constructing processes with given marginals, but none is quite suited to the problem we consider 
here. In particular, in this article we provide a family of fake exponential Brownian motions: we give 
a martingale stochastic process with continuous paths and with lognormal marginals matching those of 
exponential Brownian motion, which is not exponential Brownian motion. The ideas apply to other time- 
homogeneous diffusions. 
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One reason for interest in fake processes in this sense comes from mathematical finance, see Hobson [TT] . 
In finance knowledge of a continuum of option prices of a fixed maturity (one for each possible strike) is 
equivalent to knowledge of the marginal distribution of the underlying financial asset under the pricing 
measure. This result dates back to Breeden and Litzenberger [3J. Knowledge of the prices of options for the 
continuum of strikes and maturities is equivalent to knowledge of all the marginals of the underlying financial 
asset, but implies nothing about the finite-dimensional distributions. (However, general theory tells us that 
the asset price, suitably discounted, must be a martingale, again under the pricing measure.) In finance we 
often have knowledge of option prices, but we have little reason to believe in any particular model. The 
canonical model in finance is the Osborne-Samuelson-Black-Scholes exponential Brownian motion model, 
under which vanilla option prices are given by the Black-Scholes formula. The fact that we can provide 
alternative martingale models with the same marginals, and therefore which support the same option prices, 
shows that there is no unique model which is consistent with those prices. 

Let 5 X denote the unit mass at x. Suppose P = (Pt)t>o is exponential Brownian motion, scaled such 
that Pq = 1, and with drift parameter chosen such that P is a martingale. Setting volatility equal to one 
we may assume that P solves dPt/Pt — dB t . Let {vt)t>o be the family of laws of (Pt)t>o with lognormal 
densities (pt (•))«>() where 

1 e -*/8 e -(lnx) 2 /2t 

Ptix) = 7^^—^- (1) 



We write v t ~ C(Pt) ~ Pt{')- Following Hamza and Klebaner we ask does there exist a continuous martingale 
with marginals (i't)t>o, which is distinct from exponential Brownian motion? The main goal of this note is 
to show that the answer to this question is positive. 

For the general diffusion case, let X = (X t )t>o be a time-homogeneous martingale diffusion so that X 
solves dX t — a(X t )dB t subject to Xo = xq. Let I C R be the state space of X and let /if , with associated 
density /*(•)> be the law of X t . Assume that the boundary points of / are not attainable so that / = (£,r) 
where — oo < I < xo < r < oo and that er is positive and continuously differentiable on the interior of /. 
Further, we suppose that ft{x) is continuously differentiable in both x and t. This is automatically satisfied 
for Brownian motion and exponential Brownian motion and by results of Rogers |20) is also satisfied if a is 
thrice continuously differentiable and J e -^j^dx = oo = J r -^jdx. 

Our goal is to describe a fake version of X, where a fake version of the martingale X is a continuous 
martingale X — (Xt)t>o such that C{Xt) ~ \x£ , which is distinct from X. 

Using Jensen's inequality and the martingale property of X we can deduce that the family (fi^) t >o must 
be increasing in convex order. Further, the results of Kellerer [T3] imply that the fact that a family of 
measures is increasing in convex order is necessary and sufficient for there to be a martingale with those 
marginals. Of course, such a process need not be a diffusion. For a time-homogeneous diffusion we have the 
following result, which follows from Tanaka's formula (Revuz and Yor [l9l Theorem 6.1.2]), the fact that 
J I{x 3 >x}dX s is a martingale, and an application of Fubini's Theorem. 
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Lemma 1.1 (Carr and Jarrow [3], Klebaner [15]). Suppose er(.) is continuous. Then 

E[(X T - x)+] = (x - x) + + ^ J f t (x)dt 

2 Fake solutions based on Skorokhod embeddings 

Given a stochastic process X on a state space E, and a probability law fi on E, the Skorokhod embedding 
problem [23] is to find a stopping time t such that X T ~ ji. In the simplest Brownian case, the problem 
becomes, given Brownian motion null at zero and a centred, square-integrable probability measure fi on R, 
to find a stopping time r such that B T ~ /i and E[t] is finite (and then necessarily E[r] = J x 2 /j,(dx)). 

There are a multiplicity of solutions of the Skorokhod embedding problem for Brownian motion. Many of 
these solutions are based on considering the bivariate process (B t , A t )t>o where A is some increasing additive 
functional, null at 0. Then the solution of the embedding problem is to take r = t a — inf{i > : (B t ,A t ) e 
T> A } where 2? M is some appropriate domain in R x R + . In particular the solutions of {Root [21], Rost [22] . 
Azema-Yor [2], Vallois (24J [25]} are based on the choices A. equals {time, time, the maximum process, the 
local time} respectively. 

These constructions are suggestive of a direct construction of a discontinuous fake martingale diffusion 
X = (X t )t>o- Fix an additive functional A and let T> t = T> A X denote the domain such that r t = t a x is 
a solution of the Skorokhod embedding problem for ^ in Brownian motion. Then provided the stopping 
times r t are non-decreasing in t, or equivalently provided the stopping domains T> t are decreasing in t, the 
process X given by X t = B Tt is a martingale with the required distributions. All that remains is to check 
that the regions T> t are indeed decreasing in t. 

In general, the fact that a family of distributions (/uf r )t>o is increasing in convex order is not sufficient 
to guarantee that the associated domains T> t are decreasing in t, and for any specific additive functional and 
associated construction of the solution to the Skorokhod embedding problem this needs to be checked. 

Madan and Yor [17] use this approach and the Azema-Yor solution to describe a discontinuous fake 
diffusion. The condition that V t is decreasing in t is restated as the fact that the marginal distributions 
are increasing in residual mean life order (equivalently, the barycentre functions are increasing in t). It 
can be checked that this property holds for the lognormal family of distributions so that the Madan and 
Yor [T7] construction gives a fake exponential Brownian motion. Similarly, the Root, Rost, and Vallois 
constructions all extend from the univariate case for a single marginal, to give fake exponential Brownian 
motions. However, in all cases the resulting process is discontinuous. 

3 Continuous Fake exponential Brownian motion 

Let X — (X t )t>o be the solution of an ltd stochastic differential equation, and let nf ~ C(X t ). Then 
Gyongy [8] showed that there is a diffusion process Y = {Y t )t>o such that C(Y t ) ~ fjbf- . In a related 
result, Dupire [5] showed that if the family of distributions (nt)t>o is increasing in convex order, and if the 
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associated call price functional C(t,x) = J M (y — x)^ t {dy) is sufficiently regular (in particular, C is C 1 ' 2 and 
has certain limiting properties for large x) then there is a time-inhomogeneous martingale diffusion Y solving 
dY t = T}(t,Y t )dW t such that C{Y t ) ~ \i t . The coefficient n is given by n(t,yf = 2^^/^§^. We call 
the process Y associated with the family (nt)t>o the Dupirc diffusion. 

Given a process whose univariate marginals coincide with those of exponential Brownian motion, both 
Gydngy and Dupire show how to construct a diffusion process whose marginals are lognormal. However, 
in both cases the resulting process is exponential Brownian motion itself. Instead, in this note we give a 
family of continuous martingales which are not exponential Brownian motion but which share the lognormal 
univariate marginal distributions of exponential Brownian motion. The process we construct is a mixture of 
two well-chosen martingales. 

Let X — (X t )t>o be some time-homogeneous martingale diffusion such that X = x and for t > 0, 
C(X t ) ~ /((•)• The idea is to try to write ft(x) = cg t (x) + (I — c)h t (x) for a constant c <G (0, I ) and a pair 
of families {(gt)t>a, (ht)t>o} such that for each t, g t and h t are the densities of random variables with mean 
x , and such that if fif (respectively ^) is the law of a random variable with density g t (respectively h t ), 
and if [Iq = 5 XQ = Hq , then the family (/if ) t > (respectively (fip) t >o) is increasing in convex order. Then, 
if G = (Gt)t>o (respectively H = (Ht)t>o) is the Dupirc diffusion associated with the family of laws (fj,f)t>o 
(respectively (/i^)t>o) and if X — (X)t>o is defined by Xt — Gtl{z^=i} + Htl{z c =a} then X is a fake version 
of X . Here Z c is a Bernoulli random variable which is independent of G and H, is known at time 0, and is 
such that ¥(Z C = 1 ) = c, where c e (0, f). 

It remains to show that we can construct the families {{gt)t>o, {ht)t>o} and that the resulting process is 
not X itself. For this we need the family (/if) t >o to be increasing (in convex order), but not too quickly as 
that will mean there is no 'room' for the family (/i^) t >o to be increasing. 

The final simplifying idea is to suppose that G is a time change of X, Gt — X a ny By choosing a to be 
an increasing, differentiable process we automatically get the existence of G. If further wc require that a < 1 
we can hope that the family ()i^)t>o is also increasing in convex order. 

Proposition 3.1. Suppose a(t) is a strictly increasing, twice continuously differentiable function, null at 0, 
with d(t) < 1 for t > 0, and define 

K= inf inf ftiy \. 
Suppose K > 0. Fix c G (0, K) and for t > set 

ht{y) = -^— c {.ft{y) - cfa(t){y)} ■ 

Then (h t )t>o is a family of densities which is increasing in convex order. 
Moreover, there exists a martingale diffusion H such that C{H t ) ~ h t {-). 

Proof. From the definition of c we have that h t is non- negative and it is immediate that h t {y) integrates to 
one so that h t is the density of a continuous random variable. Further, since f t corresponds to a mean x 
random variable, so does h t . 
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Let (F f )oo be a family of random variables such that C(H t ) ~ h t (-). By Lemma [PI and the hypotheses 
of the proposition, 

^E[(H t -x)+] = J—!LE[ {Xt -x) + -c{X a{t) -x) + ] (2) 
1 - ^2 



— a{xf (f t (x) - ca(t)f a{t) (x)) > 0. 



Hence the family (/ij)t>o is a family of densities which is increasing in convex order. 
Define 

7?(*,y) =cr(y) . . t — 7^—- 3) 

JtU/) ~ cja(t)\y) 

and let iJ be a weak solution of dH t = rj(t, H t )dB t subject to Hq = xq. By our assumptions the solution to 
this SDE is unique in law. Moreover we have 

j < #,?/) 2 /t(y) < L 2 

°{y) 2 My) - c.f a (t){y) ~ 

where L 2 = K/(K — c). Further, X L given by X^ = X^i t is a martingale and solves dXf 1 — a^X^dBt 
where ctl(x) = La(x). Then, r](t,y) < (tl{v) and by Theorem 3 of Hajck [10] . we can find an increasing 
time-change T with T t < t such that H t = X^ solves dH t = rj(t, H t )dB t . Hence, by the uniqueness in law 
of solutions to this equation, H is a martingale. 

Suppose that X is non-negative (the case where X is bounded above or below reduces to this case after 
a reflection and/or a shift.) Define Ch ■ [0, oo) x [0, oo) R via Ch(0,y) = (xo — y) + and for f > 0, 
Ch(t,y) — J{z — y) + h t {z)dz. Then, using ([2} we deduce that Ch solves Dupire's equation 

1 -r 1 {t,yfC"{t, y ) = C{t,y) (4) 

subject to 

C(0, y) = (x - y)+, (7(0, t) = x , < C < x . (5) 

Indeed, by an argument similar to that in Ekstrom and Tysk [6] Step 5 of the proof of Theorem 2.2] Ch is 
the unique solution of (UJ) satisfying ([5]). 

However, by the results of Dupire [5] and Klebaner [TS] Cff(t, y) = E[(i? t — y) + ] also solves dU) and (O 
and hence, C(H t ) ~ ^-t(-)- 

Now consider the case where X the range of X is the whole real line. Let r\ and H be as before. Fix 
T e (0, oo). Let C h : [0, T] x K i-> R be the solution of gj subject to 

= (^o - y) + , sup t<T lim yt±00 C t (t,y) - (x ~ y) + = 0, (x - y) + < C h (t,y) < (xq - y)+ + J, (6) 



where J = J(T) = e ^ x l^t_ x °> By a small modification of the arguments of Ekstrom and Tysk [5] to 
allow for the fact that we are working with real- valued processes, it follows that the solution to this equation 
is unique. But Cn{t, y) = E[(i?t — y) + ] also solves (|4]) and (x — y) + < Cn(t, y) < E[(X L 2 t — y)+] and hence 
C(Ht) ~ ht for t < T. Since T is arbitrary, the result follows for all t. □ 
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Theorem 3.2. Let G be given by Gt — X a ( t ) an d let H be the Dupire diffusion associated with the family 
of densities (ht)t>o- Then 

X t = G t I{z^=i] + HtI{zc=o} 

is a fake version of X . 

Proof. If ft (x) denotes the time-i density of X we have that 

ft(x) = cg t (x) + (1 - c)ht(x) = cf a(t) {x) + (1 - c)h t {x) = f t (x). 

and Xt has the same distribution as X. 

Since G and H are martingales, so is X. But the quadratic variations of X and X are different and hence 
X is a fake version of X. □ 

Example 3.3 (Fake Brownian motion). In this case 



fa( t) (y) V t -V t ■ 

Fix c € (0, 1) and choose a{t)/t > c 2 with d(t) < 1, for example a(t) — K 2 t for some K £ (c, 1). 
Example 3.4 (Fake exponential Brownian motion). Recall the density given in {IJj. Then 



min =min ./^ e [a(t)-t]/8 e (lnx) 2 [l/a(t)-l/t]/2 = /f^) g[a(t)-t]/8 = ^( a W) 

x >° fa(t)(x) x>0\j t ' St 1p(t) 



where if) is the increasing function ip(t) = v^e*/ 8 . Fix K e (0, 1) and set a(t) = ^(Kip^)) . Then a(t) < t 

St 



and a(t) = 4>{t) / 4>{a{t)) where <fi is the decreasing function 4>{t) = ^r, so that a < 1 for t > 0. 



Remark 3.5. The argument can be extended to cover the case where the original process X is a time- 
inhomogeneous martingale diffusion dX t — a(t, X t )dB t provided k > where 

k := inf inf 1 -r—^)—- A MM*,*) 3 



*>o»ei \fa(t){y) fa{t)[y)o-{a{t),yf 
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